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ABSTRACT 


The  geomagnetic  electrokinetograph  (GEK)  was  used  to  measure 
surface  currents  near  the  center  of  Monterey  Bay  during  six  separate 
24-hour  periods  from  May  through  July,   1972.    An  average  of  244  current 
vectors  were  derived  for  each  cruise.    The  mean  currents  from  these 
cruises  are  all  southerly  and  ranged  from  4.1  cm/sec  to  20.4  cm/sec. 
The  average  of  these  mean  currents  is  12.1  cm/sec  toward  163 °T.    These 
values  were  compared  with  individual  currents  derived  from  dynamic 
topographies  from  the  same  period. 

Diurnal  and  semi-diurnal  variations  of  the  current  were  studied 
after  subjecting  the  data  to  a  fourier  analysis.     It  was  concluded  that 
there  must  be  at  least  an  indirect  coupling  of  the  ocean  currents  with 
the  semi-diurnal  tide  at  the  data  point.     The  diurnal  component  also  is 
important;  it  may  be  tidal  or  inertial,  or  merely  related  to  the  passage 
of  the  sun. 
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I.     INTRODUCTION 

I.  INTRODUCTION 

This  study  was  designed  to  use  the  Geomagnetic  Electrokinetograph 
(GEK)  to  investigate  the  time  variability  of  surface  currents  at  one  point 
in  Monterey  Bay. 

Recent  studies  by  Smith  (1971)  and  McKay  (1970)  utilized  the  GEK 
to  measure  surface  currents,  but  these  measurements  were  taken  singly, 
over  large  areas,  and  at  different  times.     In  order  to  determine  the  nature 
of  the  surface  currents  in  Monterey  Bay  it  was  necessary  to  get  data 
which  was  independent  of  spacial  changes  and  of  a  sufficiently  long 
duration  to  show  tidal  and  inertial  effects.    The  technique  followed  was 
to  continually  run  squares  with  the  GEK  for  approximately  24  hours.    The 
squares  averaged  0.5  nautical  miles  on  a  side  and  were  repeated  about 
every  24  minutes.     Electrode  zeros  and  current  components  were  com- 
puted cyclically  around  the  square. 

Six  experiments  were  conducted  between  10  May  and  22  July,    1972 
at  a  point  near  the  center  of  Monterey  Bay  (36°46'N/122  °01'W)  in  water 
440  fathoms  deep  (Figure  1).    The  Naval  Postgraduate   School  (NPS) 
research  vessel  ACANIA  was  used  to  tow  the  GEK  array.    Analysis  of 
the  data  was  facilitated  by  the  use  of  the  NPS  IBM  360  computer.     Causes 
of  variability  of  surface  currents  were  sought  in  tidal,  wind,  inertial, 
and  seasonal  effects  . 


Figure  1.  The  Location  of 
the  Data  Point  in  Monterey 
Bay. 
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II.     THE  ELECTROMAGNETIC  METHOD  OF 
MEASURING  QC EAN  C URRE NTS 


A.  DEVELOPMENT 

When  a  part  of  the  ocean  moves  relative  to  the  earth's  magnetic 
field,  a  voltage  is  induced  in  the  sea  water.    Under  certain  conditions 
a  voltage  equal  to  the  induced  voltage  may  be  measured  between  two 
identical  electrodes  introduced  into  the  sea.     Under  other  conditions, 
part  or  all  of  the  induced  voltage  is  dissipated  through  resistance  in  the 
sea  water. 

Faraday  is  credited  with  the  first  mention  of  and  experimentation 
with  this  phenomenon.    His  experiments  failed;  however,  in  the  late 
1800's  broken  trans -Atlantic  cables  showed  the  potential  differences 
he  had  sought  and  led  to  the  early  2  0th  century  experimentation  in 
Dartmouth  Harbor  by  Young,  Gerrard,  and  Jevens   (1920). 

The  electromagnetic  method  of  towed  electrodes  is  a  continuation 
of  this  theory  and  is  a  result  of  work  done  between  1946  and  1954  by 
Longuet-Higgins  and  Barber  (1946),  Stommel  (1948),  von  Arx  (1950), 
and  Longuet-Higgins,  Stern,  andStommel  (1954).    The  technique  owes 
its  name,  Geomagnetic  Electrokinetograph  (GEK)  ,  to  von  Arx.     Using  the 
GEK  method,   surface  currents  can  readily  be  measured  from  a  ship  underway, 

B.  GEOMAGNETIC  ELECTROKINETOGRAPH   (GEK) 
1 .         Theory  of  Operation 

The  GEK  is  made  up  of  a  recording  potentiometer,  a  two- 
conductor,  neutrally  buoyant  cable,  and  a  matched  pair  of  silver-silver 
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chloride  electrodes  of  the  order  of  100  meters  apart.    The  potentiometer 
records  the  difference  in  potential  between  two  electrodes.    A  long  cable 
is  required  to  allow  the  electrodes  to  be  towed  beyond  the  magnetic 
effects  of  the  ship's  steel  hull. 

An  extensive  discussion  of  electromagnetic  theory  and  the  ocean's 
movement  through  the  earth's  magnetic  field  may  be  found  in  von  Arx 
(1950)  and  Longuet-Higgins ,  Stern,  and  Stommel  (1954)   (hereafter  referred 
to  as  L-HSS)  .     The  following  summary  of  L-HSS's  treatment  of  this  sub- 
ject is  by  McKay  (1970)  with  large  portions  of  the  development  based  on 
work  by  Curtin  (1970). 

The  electrical  field  intensity  along  a  conducting  path  that  moves 
through  a  magnetic  field  is  the  difference  between  the  two  potential 
gradients.     One  potential  gradient  can  be  shown  by  application  of 
Faraday's  Law  of  Induction  to  be  V  x  H  as  is  done  by  L-HSS.     In  this 
expression,  V  is  the  vector  velocity  with  which  the  conducting  path  is 
moved  and  H  is  the  magnetic  field  intensity  vector.     The  other  potential 
gradient  is  given  by  Ohm's  Law  as         P ±->  where        a    is  the  resistivity 
and_i_the  electric  current  density  vector  along  the  conductor.     If  the 
potential  is  denoted       <E>       ,  the  negative  of  the  electric  force  field  in 
the  moving  conducting  path  can  be  expressed  as 

V$  =  (VxH)  -        >0i  (1) 

When  the  electrodes  are  deployed,  two  distinct  conducting  paths 
are  present:  (a)  the  sea  water,  and  (b)  the  long  conductors  to  which  the 
electrodes  are  attached.    The  potential  in  each  of  these  paths  is  the 
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difference  of  the  two  terms  of  equation  (1) .     In  the  following  discussion, 
a  right-hand  co-ordinate  system  is  utilized  with  the  following  conven- 
tions:   the  positive  x-direction  is  that  of  the  ship's  head;  y  is  the 
horizontal  co-ordinate  perpendicular  to  the  course  steered  and  z  is 
positive  up. 

Consider  first  the  sea  water  path.    When  the  electrodes  are  towed 
directly  astern  of  the  ship,  as  shown  in  Figure  2,  they  sense  the  poten- 
tial difference  along  a  path  parallel  to  the  direction  of  the  ship's  head. 
The  component  of  potential  gradient  generated  in  this  direction  due  to 

the  (V  x  H)  term  is   (V  H     -  V  H   ),  where  V     is  the  component  of  water 

y    z  z    y  y 

velocity  perpendicular  to  the  course  steered  (set  of  ship),  V     is  the 

z 

vertical  component  of  water  velocity,  H     is  the  vertical  intensity  of  the 

z 

earth's  magnetic  field,  and  H  is  the  horizontal  intensity  of  the  geo- 
magnetic field  in  a  direction  perpendicular  to  the  ship's  course.     If  V 

z 

is  assumed  negligible  in  comparison  with  V    ,  the  measurable  potential 

y 

gradient  generated  by  the  water  motion  is  V  H    ;  that  is,  the  potential 

y   z 

gradient  contributed  by  the  (VxH)  term  is  due  to  the  component  of  water 
velocity  perpendicular  to  the  ship's  course.    This  water  motion  causes 
the  ship's  set  (and  the  set  of  the  electrodes)  and  is  independent  of  the 
ship's  velocity  in  the  direction  of  the  interelectrode  line.     The  second 
potential  gradient  in  the  sea  water  path,  /O  J_,  is  a  direct  function 

of  the  electric  current  density  along  the  path  between  the  electrodes. 
The  electric  current  density  along  this  path  depends  on  the  electrical 
characteristics   (potential  gradient's  resistances)  of  all  other  surrounding 
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paths,   since  it  is  only  through  the  surrounding  paths  that  a  closed 

electric-current -carrying  circuit  can  be  completed. 

Consider  next  the  part  of  the  closed  conducting  path  formed  by  the 

two  wires  and  the  voltmeter.    It  can  be  seen,  by  application  of  reasoning 

identical  to  that  used  for  the  sea  water  path,  that  the  velocity -induced 

potential  gradient  due  to  the  (V  x  H)  term  in  equation  (1)  is  again  V  H    . 

-  y    z 

The  second  potential  gradient        fi ±.>  is  essentially  zero  in  this  path 
due  to  the  low  resistance  of  the  wire  and  the  small  current  drawn. 

The  potential  differences  sensed  by  the  voltmeter  and  continuously 
recorded  as  GEK  data  are  those  developed  in  the  wires  leading  to  the 
electrodes  and  in  the  sea  water  between  the  electrodes.    Traversing  this 
path,  as  shown  in  Figure  2,  and  adding  the  potentials  yields  the  following: 


E  =  (V  H   )AB  +  (V  H   )BC  +      [(-V  H   )CB  -     (    />  i   )CB] 
yz  yz  y    z  /       n 

+  (-V   H   )BA 
y    z 


(2) 


where  E  is  the  emf  per  unit  distance  and  A,   B,  and  C  are  the  points  shown 
in  Figure  2.    The  sumbol  AB  represents  the  distance  between  A  and  B. 
The  term  in  brackets  represents  the  potential  gradient  in  the  sea  water 
between  the  electrodes;    the  remaining  terms  represent  the  induced  emf 
in  the  indicated  sections  of  the  wire.    Equation  (2)  reduces  to: 

E  =  (      jO  i   )CB  (3) 

where  the  direction  is  that  of  the  electrode  line,  or,   in  the  case  of 
electrodes  towed  directly  astern,  the  direction  of  the  ship's  head. 
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For  most  currents  in  the  deep  ocean  the  potential  gradient  is 
completely  short-circuited  by  the  deep  layer  of  relatively  still  water 
underlying  the  moving  surface  layer.    In  this  case         \7  CD    -  0  and 

/2-i  ~  Y.x  S-     Under  these  conditions,  if  the  interelectrode  distance 
is  s  units,  the  horizontal  potential  difference  measured  by  the  voltmeter 
is: 

E  =  (ksHzVh)  xlO    "  (4) 

where  V,    is  the  horizontal  component  of  water  motion  perpendicular  to 
n 

the  interelectrode  line  and  k  is  a  constant  near  unity,  to  be  discussed 

in  Section  III  C . 

The  positive  direction  of  E  is  taken  as  that  of  the  ship's  head 

and  the  current  velocity  is  determined  by  adding  the  electric  vectors,  E, 

obtained  on  two  courses,  at  right  angles.    The  following  example  indicates 

the  relationship  between  the  direction  of  the  current  component  measured 

and  that  of  the  interelectrode  line: 

If  V     is  positive  toward  the  east,  H     is  directed  downward, 
h  z 

and  the  interelectrode  line  is  parallel  with  a  meridian  of 

longitude,  E  will  have  a  positive  polarity  at  the  southern 

electrode  and  a  negative  polarity  at  the  northern  electrode. 

Since  H     is  positive  downward  in  the  northern  hemisphere,  the 
z 

following  convention  must  be  observed  to  obtain  the  true  direction  of 
the  current  velocity:    In  order  that  the  deflection  recorded  by  the  in- 
strument be  in  the  same  direction  as  the  set  of  the  electrodes,  the 
conductor  leading  to  the  more  distant  electrode  must  be  connected  to 
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the  voltmeter  input  terminal  which,  when  made  positive,  gives  a  de- 
flection of  the  recorder  current  indicator  to  the  right.    The  necessity  of 
observing  this  convention  can  be  seen  from  the  example  given  above. 
This  convention  must  be  reversed  in  the  southern  magnetic  hemisphere, 

where  H     is  directed  upward. 
z 

2  .         Assumptions 

As  has  been  stated  above,  the  electric  current  density,  j_, 
is  a  function  of  the  conducting  path  in  which  it  exists .    Thus ,  with  an 
assumption  of  uniform  resistivity,  the  current  density  integrated  along 
any  line  is  directly  dependent  on  the  emf's  in  the  available  sea  water 
paths  surrounding  that  line.     These  emf's  as  indicated  by  Equation  (1), 
are  directly  related  to  the  water  velocity  in  the  region.    Therefore,  any 
electric  current  density  measured  at  the  surface  is  a  reflection,  in  part, 
of  the  velocity  field  in  the  entire  water  column  below  the  surface  and  in 
adjacent  waters  . 

Recognition  of  this  fact  is  essential  for  a  meaningful  inter- 
pretation of  GEK  data.     It  is  only  in  the  case  of  true  oceanic  currents 
that  GEK  data,  unsupplemented  by  information  derived  from  other  sources, 
can  be  considered  representative  of  the  true  surface  current  velocity.     In 
this  case  the  following  conditions  generally  prevail: 

(i)    the  vertical  component  of  water  velocity  is  negligible 
in  comparison  with  the  horizontal  component, 

(ii)    the  flow  is  broad  compared  with  the  depth  of  the  water, 

(iii)    the  current  is  steady  in  the  x-diroction,  while  a  current 

"fix"  is  being  obtained, 
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(iv)    a  deep  layer  of  still  water  underlies  the  moving  surface 
layer. 
If  condition  (iv)  is  not  obtained,  the  voltage  is  usually  less  than 
the  theoretical  and  it  is  customary  to  multiply  by  a  correction  factor  k 
determined  by  a  simultaneous  measurement  with  a  current  meter. 
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III.     EQUIPMENT 

A.  ELECTRODES 

A  complete  operating  assemblage  of  GEK  equipment  was  inherited 
from  Smith  (1971)  at  the  onset  of  this  project. 

The  matched  set  of  silver-silver  chloride  electrodes  was  supplied 
by  Woods  Hole  Oceanographic  Institute  and  attached  to  the  cable  100 
meters  apart  by  direct  soldering.    Von  Arx  (1950)  describes  the  process 
by  which  electrode  pairs  are  matched.    The  electrodes  were  padded  and 
housed  in  a  water  permeable  tube  as  shown  by  Smith  (1971). 

During  the  three-month  period  of  data  collection,  it  was  not  neces- 
sary to  change  the  electrode  installation. 

To  eliminate  the  electrochemical  potential,  the  "zero  point"  must 
be  obtained  and  subtracted  out  of  the  actual  current  vectors.     This  is 
normally  done  by  reversing  the  electrodes  end  for  end.     In  this  case,  the 
ship  used  a  square  sailing  plan  and  the  zero  point  was  taken  from  opposite 
sides  of  each  square.    This  is  described  in  detail  in  Section  V  A. 

B .  CABLE 

A  neutrally  buoyant  cable  (Type  B-3088)  329  meters  long,  made  by 
the  Vector  Cable  Company,  was  utilized  during  this  investigation.    The 
cable's  function  is  to  transmit  to  the  recorder  the  electrode-generated 
signals  from  far  enough  astern  (1.5  -  3.0  shiplengths)  of  the  ship  that 
the  electrodes  are  not  affected  by  the  ship's  magnetic  field  and  wake. 
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C.         RECORDING  POTENTIOMETER 

The  Geomagnetic  Electrokinetograph  Model  V  recording  potentiometer 
based  on  a  Speedomax  Type  G  recorder  was  used  to  record  the  potential 
differences  delivered  from  the  electrodes  by  the  cable.    A  photograph 
and  a  thorough  description  of  the  operation  of  the  GEK  Model  V  is  con- 
tained in  H .  O.  Publication  60  7. 

The  vertical  intensity  of  the  earth's  magnetic  field  in  gauss  for 
the  ship's  location  was  taken  from  H.  O.  Chart  Number  1702  and  the 
rounded  value  was  set  into  the  recorder  to  the  nearest  tenth.    A  multi- 
plication factor  was  later  used  to  adjust  the  data  to  the  exact  magnetic 
vertical  intensity  at  the  data  point.    The  multiplier  includes  compensation 
for  annual  variation,  and  is  the  ratio  of  the  rounded  vertical  intensity  to 
the  correct  local  magnetic  intensity.    In  this  instance  it  was  0.91 
(0.40/0.44). 

As  noted  by  Smith  (1971)  it  is  essential  to  isolate  the  power  for 
the  recorder  from  the  ship's  ground.    This  is  done  with  an  isolation 
transformer  which  must  not  be  by-passed  by  the  ground  lead  of  a  3-wire 
power  cord.     Hence,  a  conventional  two-prong  plug  is  necessary. 

The  recorder  also  requires  a  signal  multiplier  setting.    This  k 
value  is  the  ratio  of  the  water  speed  observed  by  a  non-electromagnetic 
method  to  the  water  speed  measured  by  the  electromagnetic  method.    An 
open  ocean  value  of  1.04  is  suggested  by  von  Arx  (1950)  with  the  caution 
that  higher  values  might  be  needed  in  areas  where  tidal  currents  and 
water  shear  occur. 
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The  k  factor  was  set  at  1 .  10 ,.  following  the  recommendation  of 
H.  O.  Publication  607  (1970  Ch.   1)  for  water  deeper  than  75  fathoms. 
This  factor  was  derived  from  experience  in  open  ocean  conditions  along 
the  Atlantic  coast.    An  unusual  structure  like  Monterey  Canyon  may 
behave  anomalously,  a  possibility  which  will  be  considered  in  the  dis- 
cussion of  accuracy  (Section  V  D) .    An  incorrect  choice  of  k  does  not 
lead  to  an  error  in  direction. 
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IV.     SHIP  PROCEDURES 

A.  PLATFORM 

The  research  vessel  ACANIA  was  acquired  by  the  Naval  Postgraduate 
School  in  1971  and  refitted  as  an  oceanographic  research  vessel.     She 
was  launched  in  1929  and  for  many  years  took  part  in  research  activities 
under  the  Stanford  Research  Foundation. 

R/V  ACANIA's  length  is  126  feet,  maximum  speed  is  ten  knots,  and 
she  carries  a  crew  of  five.     Of  the  standard  equipment  aboard  ACANIA, 
the  Decca  RM  416  radar  was  the  most  important  to  this  study  as  noted  in 
Section  IV  C.    The  UQN  -  1  depth  finder  was  also  used  briefly  to  confirm 
charted  depths  at  the  data  point. 

B.  SAILING  PLAN 

Both  clockwise  and  counterclockwise  squares  were  traversed  during 
the  cruise  series,  allowing  the  ship  to  remain  within  1000  meters  of  the 
data  point  in  calm  weather  and  within  1500  meters  during  heavier  weather. 

Ship's  speed  was  held  at  five  knots  for  all  cruises,  but  the  time 
on  each  leg  was  varied  between  four  and  nine  minutes,  depending  upon 
wind  and  current  conditions.     In  order  to  maintain  position  within  1000 
meters  of  the  data  point  the  downwind  leg  might  take  four  minutes,  while 
the  upwind  leg  (opposite  side  of  the  square)  might  take  eight  minutes. 
All  legs  were  at  least  four  minutes  long  to  insure  that  the  cable  had 
straightened  out  and  a  good  signal  was  being  received. 
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FIGURE  3.     .SAILING  PLAN  AND  COMPONENT  COMPUTATION 


FIRST  VECTOR: 


ZERO  =  VI  +  V3 


E-W  COMPONENT  =  Vl-ZERO 
N-S  COMPONENT  =  V2-ZERO 

2  2     1/2 

CURRENT  MAGNITUDE  -  [  (E-W)     +  (N-S)     ] 

CURRENT  DIRECTION  =  arcsin  E-W 

MAGN 


SECOND  VECTOR: 


ZERO  =  V2  +  V4 


E-W  COMPONENT  =  V3-ZERO 
N-S  COMPONENT  =  V2-ZERO 


THIRD  VECTOR: 


ZERO  =  V3  +  V5 


E-W  COMPONENT  =  V3-ZERO 
N-S  COMPONENT  =  V4-ZERO 
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Only  cardinal  headings  were  steered  during  this  series  of  cruises. 
The  sailing  plan  is  shown  in  Figure  3  . 

C .         NAVIGATION 

The  data  point  for  this  project  was  located  close  to  the  center  of 
Monterey  Bay  over  the  Canyon.     Navigation  was  by  radar  at  ranges  of 
roughly  11  miles  from  well-defined  landmarks.     These  points  are  noted 
on  Figure  1.    Visual  bearings  were  also  used  to  augment  the  radar  range 
information.    The  positioning  accuracy  is  estimated  to  be  +  0 . 1  nm. 
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V.     DATA  COLLECTION,  REDUCTION,  AND  PRESENTATION 

A.         SURFACE  CURRENTS 
1 .         Data  Collection 

During  each  cruise  approximately  twenty-four  hours  of  con- 
secutive surface  current  readings  were  obtained.    All  data  was  taken  at 
one  point  near  the  center  of  both  Monterey  Bay  and  the  Monterey  Canyon. 
An  average  of  244  current  vectors  were  obtained  from  each  of  the  six 
cruises . 

The  GEK  recording  potentiometer  was  continuously  manned 
in  order  that  the  strip  chart  could  be  annotated  and  the  potentials  logged 
for  each  leg  of  the  sailing  plan.     The  ship's  speed  was  five  knots  with 
each  leg  of  the  square  being  between  four  and  nine  minutes  in  duration. 
The  average  duration  was  six  minutes.     Longer  and  shorter  legs  were 
run  to  compensate  for  the  ship's  drift. 

The  primary  items  logged  on  the  data  sheet  were:    ship's 
heading  (°T) ,  velocity  (cm/sec),  and  time.    Weather,  position,  depth, 
speed,  and  other  pertinent  data  were  also  logged  on  each  data  sheet. 

The  strip  chart  was  marked  and  the  time  was  noted  for  each 
change  of  course.  The  date,  segment  of  cruise  and  cruise  number  were 
also  recorded  on  the  chart. 

Breaks  in  the  data  of  about  ten  to  fifteen  minutes  occurred 
when  the  direction  sailed  around  the  square  was  changed  from  clockwise 
to  counterclockwise  or  vice  versa.    Several  breaks  in  the  data  sets  of 

up  to  forty-five  minutes  occurred  due  to  steering  casualties. 
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2  .         Data  Reduction 

The  IBM  360/67  computer  at  the  Naval  Postgraduate  School 
computer  center  was  used  to  reduce  the  rough  data  to  current  magnitudes 
and  directions . 

As  can  be  seen  from  Figure  3  ,  there  is  an  interdependence 
of  current  vectors  in  each  square  since  one  of  the  current  vector  com- 
ponents is  calculated  from  a  leg  of  the  square  which  will  take  part  in  the 
next  zero  computation.     This  procedure  allowed  for  the  calculation  of 
approximately  ten  current  vectors  per  hour. 

The  outputs  of  the  data  reduction  computer  program  included 
east-west  (E-W)  and  north-south  (N-S)  components,  time,  and  the  mag- 
nitude and  direction  of  each  current  vector.    This  computer  program  is 
listed  in  Appendix  A  with  suitable  comments  and  explanations.    Appendix 
B  lists  sample  current  data  for  Cruise  I. 

Sharply  alternating  speeds  were  present  in  the  data  after 
reduction.     These  variations  were  removed  by  averaging  alternate  pairs 
of  data.    A  full  discussion  of  these  variations  is  to  be  found  in  Section  V 
D. 

3  .         Data  Presentation 

To  produce  a  clear  presentation  of  current  magnitudes  versus 
winds  and  tides,  a  computer  program  was  written  which  utilized  the 
plotting  package  for  the  CALCOMP  Model  765  plotting  equipment  along 
with  the  IBM  360  computer. 
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For  each  cruise  this  graphing  program  produced  the  plots  of  E-W 
and  N-S  components  and  resultants  versus  time.     Corresponding  tide 
heights  and  wind  speeds  were  plotted  on  the  same  graph.     Northerly  and 
easterly  sets  were  taken  to  be  positive.     Both  unaveraged  and  averaged 
data  were  handled  in  this  fashion.    The  graphing  program  is  included 
after  the  reduction  program  in  Appendix  A. 

The  data  is  also  presented  in  a  progressive  vector  diagram  format 
for  each  cruise.  The  progressive  vector  diagram  again  utilized  the  plotting 
package  and  may  be  found  with  explanations  in  Appendix  A. 

B.  WINDS 

Estimations  of  wind  speed  and  direction  were  recorded  aboard  R/V 
ACANIA  at  various  times  (normally  every  two  or  three  hours)  during  each 
cruise.  These  estimates  were  compared  with  the  hourly  readings  logged 
at  the  U.S.  Coast  Guard  lighthouse  on  Point  Pinos .  A  composite  of 
these  two  sources  forms  the  basis  for  the  wind  magnitude  and  direction 
utilized  for  each  cruise.  Appendix  C  contains  the  resulting  hourly  wind 
data  for  Cruise  I. 

C.  TIDES 

Tide  records  from  Pier  Number  1  at  Monterey  were  made  available 

by  the  NPS  Oceanographic  Department.     This  data  was  used  to  plot 

directly  the  tidal  curve  for  each  cruise.    These  tides  compare  favorably 

with  the  high  and  low  tidal  heights  as  calculated  from  the  U.  S.  Coast 

and  Geodetic  Survey  Tide  Tables.    The  hourly  tidal  heights  utilized  are 

contained  in  Appendix  D. 
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A  tide  prediction  computer  program  was  also  written,  based  on 
harmonic  constants  for  San  Francisco  and  corrected  to  the  data  point  in 
Monterey  Bay.     The  results  from  this  effort  did  not  correlate  well  with 
either  the  tide  gauge  measurements  or  the  tide  tables,  and  were,  there- 
fore, not  considered  further. 

D.         DISCUSSION  OF  ACCURACY 

In  discussing  the  accuracy  of  the  current  data  the  following  poten- 
tial errors  must  be  considered: 

(i)    Errors  in  the  observed  voltage  due  either  to  reading  in- 
accuracies or  to  the  failure  of  the  trace  to  steady  to  a 
constant  value  after  a  turn, 
(ii)    Errors  due  to  inadequate  depth  or  effects  of  the  canyon 

walls  not  compensated  for  by  the  k  factor, 
(iii)    Errors  due  to  substantial  velocity  gradients  across  the 
square . 
Since  the  trace  could  be  read  to  +  1  cm/sec,  and  at  least  four 
minutes  were  always  given  to  steady  out  this  trace,  the  errors  in  Group 
(i)  are  small.     They  are  not  systematic  and  in  the  long  term  average  to 
zero . 

It  is  not  known  how  great  the  errors  of  Group  (ii)  may  be.    There 
may  not  be  a  thick  layer  of  static  water  near  the  bottom  (current  shears, 
internal  waves  may  exist)  and  the  lateral  dimensions  of  the  parcel  of 
deep  water  are  relatively  small.     Because  of  the  sharply  sloping  walls 
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of  the  canyon,  bottom  effects  may  affect  different  legs  unequally  and 
cause  aberrations.     In  these  circumstances  it  is  unwise  to  be  categorical 
about  the  adequacy  of  using  k  =  1.10.     On  the  other  hand,  the  speeds 
measured  are  plausible  and  not  notably  at  variance  with  the  few  available 
measurements  obtained  in  various  locations  in  Monterey  Bay. 

The  sharply  alternating  speeds  described  in  Section  V  A  2  may  best 
be  explained  by  the  existence  of  large  velocity  gradients  across  the 
square.     These  gradients  may  be  caused  by  variations  of  depth  or  by 
actual  currents  or  eddies.    This  is  contrary  to  the  assumption  that  the 
current  is  broad  and  therefore  the  same  on  both  sides  of  the  sailing  plan. 
This  breakdown,  coupled  with  the  circular  and  redundant  method  by  which 
the  currents  were  computed,  produces  oscillations  in  a  sequence  high- 
medium-low-medium-high  in  the  N-S  and  E-W  components  of  current 
and  in  the  resultant.   (See  Appendix  C  for  an  example.)    If  the  current 
component  actually  present  on  one  N-S  leg  has  a  magnitude  E  and  on 
the  opposite  side  has  the  magnitude  E  +       /£\  E,  and  likewise  the  re- 
maining two  sides  of  the  square  have  magnitudes  N  and  N  +        /\  N  , 
then  it  may  be  shown  that  the  peak-to-peak  range  of  the  resulting  zig- 
zag will  be  the  absolute  value  of        A  E  -        A  N.     (Sec  Tigure  3  for 
method  of  computation.)       A  maximum  peak-to-peak  range  of  60  cm/sec 
occurs  in  the  E-W  component  of  Cruise  I,  though  most  peak-to-peak 
variations  are  of  the  order  of  30  cm/sec.    These  variations  were  mostly 
removed  from  the  data  by  averaging  alternate  pairs  of  values.    Variations 
after  averaging  are  of  the  order  of  5  cm/sec.    The  process  of  averaging 
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alternate  pairs  would  produce  a  speed  corresponding  to  some  point  within 
the  square,  if  the  first  explanation  of  the  alternation  is  valid.     It  can  be 
seen  from  comparing  the  graphs  of  averaged  and  unaveraged  data,  in 
Section  VI  and  Appendix  C,  that  the  averaging  process  produces  a  plausible, 
smoothly  varying  curve. 
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VI.     ANALYSIS  OF  CURRENT  DATA 

A.         INTRODUCTION 

Various  presentations  were  used  to  show  different  aspects  of  the 
current  data,  and  to  correlate  this  data  with  winds,  tides,  and  ocean 
currents.    A  summary  of  results  of  the  fourier  analysis  with  mean  current 
values  and  a  table  of  contemporary  current  data  are  included  in  this 
chapter  (Tables  I  and  II) . 

The  graphical  presentations  are  contained  in  Section  B  of  this 
chapter  in  the  following  order: 

(i)         East-west  current  component  versus  wind  and  tide, 
(ii)         North-south  current  component  versus  wind  and  tide, 
(iii)         Total  current  magnitude  versus  wind  and  tide, 
(iv)         Progressive  vector  diagram  with  the. cumulative  trans- 
port past  the  data  point  marked  every  five  "fixes" 
(approximately  30  minutes)  with  a  cross. 
Each  set  of  graphs  is  prefaced  with  a  short  summary  of  sea  and 
weather  conditions  for  that  particular  period. 
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B.         GRAPHICAL  PRESENTATION  OF  DATA  BY  CRUISE 
1 .         Cruise  I 

Period:     10-11  May  1972 
Latitude:     36°46'  N 
Longitude:     122° 01'  W 
Depth:    440  fathoms 
Ship's  Speed:     5  knots 

Sea  and  Weather:    Winds  were  from  the  WNW  during  this 
cruise,  with  a  diurnal  variation  in  wind  speed  caused  partly 
by  a  land-sea  breeze  pattern.    Velocities  of  up  to  22  knots 
occurred  in  the  late  afternoon  and  virtual  cairn  prevailed  in 
the  early  morning  hours.     Seas  varied  from  0.5  to  2.0  feet 
with  the  wind  velocity,  while  the  swell  was  between  3  feet 
in  the  early  morning  and  8  feet  in  the  early  evening.    Both 
were  from  the  WNW.    The  sky  was  partly  cloudy.     Note 
that  only  the  averaged  current  data  is  presented  in  this 
section.    The  unaveraged  current  data  is  presented  in 
Appendix  C.    The  start  of  the  progressive  vector  diagram 
is  labelled  "S"  and  the  end  is  labelled  "E"  . 
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Figure  5d.   Cruise  I 
Progressive  Vector  Diagram, 
N-S  Component  versus  E-W 
Component;  Averaged  Current  Data 
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Cruise  II 

Period:    30  -  31  May  1972 

Sea  and  Weather:    Same  pattern  and  magnitude  as  in  Cruise  I, 

but  with  fog  and  overcast  skies  in  the  morning,  clearing  by 

mid-afternoon. 

Note:    All  other  information  is  the  same  as  for  Cruise  I. 
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Figxu-a  6d„   Cruise  II 
Progressive  Vector  Diagram, 
N-S  Component  versus  E-W 
Component;  Averaged  Current  Data 
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3  .         Cruise  HI 

Period:     19-  20  June  1972 

Sea  and  Weather:    Weak  and  variable  winds  with  a  maximum 
velocity  of  8  knots  were  the  rule,  with  a  1  foot  sea  and  3 
foot  swell  throughout  the  cruise.     The  sky  was  overcast  and 
some  drizzle  fell  in  the  morning . 
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Figure  7d.   Cruiso  III 
Progressive  Vector  Diagram, 
N-S  Component  versus  E-W 
Component;  Averaged  Current  Data 
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Cruise  IV 

Period:     22  -  23  June  1972 

Sea  and  Weather:    Very  similar  to  Cruise  I  with  a  maximum 

wind  of  16  knots  and  overcast  skies.    Winds  varied  from 

west  to  southwest  during  this  period . 
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Figur®  8d.   Cruise  IV 
Progressive  Voctor  Diagram, 
N-3  Component  versus  E-W 
Component;  Av©rag©d  Current  Data 
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5  .         Cruise  V 


Period:     10-11  July  1972 

Sea  and  Weather:    Wind  was  fairly  steady  from  the  west  to 
northwest  at  5  to  18  knots  with  the  same  diurnal  pattern  as 
Cruise  I.    The  skies  were  clear.     Seas  and  swell  were  from 
the  WNW  with  the  same  magnitude  as  Cruise  I. 
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Figure-  9d.  Cruisa  V 
Progressiva  Vector  Diagram, 
N-S  Component  versus  E-W 
Component;  Averaged  Current  Data 
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Cruise  VI 

Period:    20-21  July  1972 

Sea  and  Weather:    Winds  were  predominantly  southerly  with 

velocities  below  10  knots.    There  was  a  mixed  swell,  but 

it  was  mostly  from  the  west  and  southwest  with  a  southerly 

sea  of  1  to  2  feet.    It  was  overcast  during  much  of  this 

cruise . 
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Flgura  lOd.   Cruise- VI 
Progressive  Voctor  Diagram, 
N-S  Component  vorsus  E-W 
Component;  Averaged  Current  Data 
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C.        DISCUSSION 
1 .         General 

This  section  deals  with  the  mean  currents  and  their  relation 
to  the  wind  and  seasonal  effects.    The  currents  in  Monterey  Bay  are 
weak  and  variable  and  have  as  yet  defied  accurate  description.     The 
current  data  gathered  during  this  research  and  its  relationship  to  oceanic 
currents,  the  winds,  and  seasonal  effects  may,  together  with  other 
investigations,  lead  to  a  viable  picture  of  the  bay's  circulation. 

The  averages  of  the  mean  components  for  the  six  cruises 
are  11.6  cm/sec  to  the  south  and  3.6  cm/sec  to  the  east.    These  yield 
a  resultant  of  12.1  cm/sec  toward  163°T.     No  data  are  available  which 
may  be  compared  accurately  to  the  present  results.    However,  for  purposes 
of  rough  comparison,  Table  I  gives  the  mean  for  each  cruise  and  currents 
derived  from  dynamic  heights  by  Moomy  (results  yet  to  be  published) . 
The  points  between  which  Moomy' s  transports  are  computed  are  shown 
in  Figure  4.     The  mean  of  Moomy' s  data  is  15  cm/sec,  and  of  the  present 
data,   12  cm/sec.    However,  the  standard  deviations  are  relatively  large 
and  the  close  agreement  probably  is  fortuitous.     However,  the  generally 
southerly  trend  of  all  the  mean  directions  of  the  time  series  data  is  in 
good  agreement  with  Moomy' s  data. 

The  "upwelling  period,"  during  which  this  current  data  was 
collected,  has  been  described  by  Skogsberg  (1936)  as  extending  from 
February  to  September,  and  as  having  the  same  prevailing  northwesterly 
winds  as  were  experienced  during  this  research.     Lammers   (1971)  noted 
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TABLE  1 


COMPARISON  OF  RECENT  OCEAN  CURRENT  MEASUREMENTS 


HOV/TON'S 

DATA 

MOO  MY' S  DATA 

MEAN  SURFACE  CU! 

[RE 

NTS  FROM 

24-HOUR  TIME  SERIES  OF 

TRANS 

.PORT  DERIVED  FROM 

GEK 

MEASUREMENTS 

DYNAMIC  HEIGHTS 

(AVERAGE  OF 

6  MEANS: 

12 

1/163°T) 

VEL/D1RC 

VEL/D1RC 

'  DATE 

STATION 
Data 

CM/SEC/°T 

DATE 
9  May 

STATION 

3  &  5 
3  &  4 

CM/SEC/"  T 

21  /  SE 
26  /SE 

10/11  May 

Point 
(Near 

7.2  /  164 

16  May 

4  &  5 
2  &  4 

2  &  3 
2  &  3 

3  /  NE 
8  /S 
27  /  S 
13  /N 

30/31  May 

Station  #5)    17.1  /  175 

17  May 

2  &  3 

2  /  S 

19/20  June 

Data 

4.1/198 

2  June 

253 

11  /S 

22/23  June 

Point 

11.3  /  135 

10/11  July 

Data 

21.9  /  154 

6  July 

2  &  3 

22  /  S 

Point 

19  July 

2  &  3 

7/S 

20/21  July 

11.9  /  158 

25  July 

2  &  3 

14  /  S 

Figure  ^.  Loc&tiono  of 
Other  Recent  Current 
Measurements* 
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that  the  southern  flow  of  the  California  Current  exists  at  the  seaward 
edge  of  Monterey  Bay  from  mid -February  to  November,  while  the  northerly 
Davidson  Current  persists  during  the  balance  of  the  year.     Lammers  also 
observed  that  during  the  upwelling  period  the  mean  velocity  of  the  Cali- 
fornia Current  in  this  region  is  from  1.0  to  6.7  nm/day  with  a  mean  of 
5  nm/day.    This  generally  agrees  in  both  direction  and  speed  with  the 
average  of  the  mean  currents  from  each  cruise  (5.9  nm/day). 

During  this  set  of  cruises  the  wind  blew  fairly  steadily  from 
the  west  northwest  with  the  exception  of  Cruise  VI,  when  the  wind  was 
from  the  south  and  southwest  at  ten  knots  or  less.     If  the  general 
southerly  flow  is  due  to  the  prevailing  wind,  this  cruise  would  be  ex- 
pected to  show  northerly  flows.    However,  during  Cruise  VI  the  mean 
current  maintained  its  mean  southerly  flow  with  only  slightly  more  than 
the  usual  easterly  component.     In  checking  the  progressive  vector 
diagram  for  Cruise  VI  it  was  seen  that  almost  due  easterly  current  flow 
occurred  during  the  daylight  hours  when  the  southerly  wind  was  at  its 
peak.    A  further  study  of  weather  charts  for  the  period  showed  that  the 
south  wind  was  only  a  local  phenomenon  and  that  west  northwest  winds 
prevailed  along  the  coast  to  the  north.    The  local  south  wind  did  have 
some  temporary  effect  on  the  surface  currents,  but  the  California 
Current  maintained  by  the  prevailing  northwest  winds  along  the  coast 
persisted  in  its  steady  southerly  flow  across  the  mouth  of  the  bay. 
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2.         Variability 

By  inspection  of  the  unaveraged  Total  Current  Magnitude 
Graphs  (Appendix  C)  it  can  be  seen  that  the  range  of  variability  over 
24  hours  is  from  2  cm/sec  to  around  40  cm/sec.    The  Cruise  V  Total 
Current  Magnitude  graph  (Figure  15c)  is  a  good  example  of  this  variation 
and  also  shows  two  unusually  high  magnitudes  of  over  50  cm/sec.    The 
spacings  of  the  crosses  on  the  Progressive  Vector  Diagram  also  give  an 
indication  of  this  current  variability  (each  cross  represents  five  current 
vectors,  or  approximately  30  minutes  of  time)  as  well  as  the  changeability 
of  the  current  direction.     The  Cruise  I  Progressive  Vector  Diagram 
(Figure  5d)  shows  the  clockwise  rotation  (which  may  vary  through  more 
than  180  degrees)  that  is  typical  of  all  cruises  except  Cruise  V.     This 
directional  oscillation  is  superimposed  on  a  moderate  mean  southerly 
flow. 

In  order  to  produce  numerical  estimates  of  the  cyclic  com- 
ponents of  the  currents  and  to  determine  their  origins  ,  the  data  was 
subjected  to  a  fourier  analysis  .     It  was  recognized  that  with  so  short 
a  series  (24.8  hours)  the  diurnal  component  would  not  be  separated  from 
the  inertial  component  (20  hours  at  this  latitude),  and  that  if  there  were 
a  substantial  amount  of  noise  it  would  have  its  usual  effect  in  reducing 
the  reliability  of  the  coefficients.     The  N-S  and  E-W  components 
separately  were  prepared  by  filling  gaps,  extrapolating  to  24.8  hours 
total  length,  and  by  subtracting  the  mean.     Fourier  coefficients,   mag- 
nitudes, and  phase  angles  were  then  computed,  using  the  Subroutine 
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FORIT  from  the  IBM  Scientific  Subroutines  Package  (included  in  Appendix 

A).    This  program  expresses  the  coefficients  A     and  B     as  those  in  the 

m  m 


equation 


f(t)=A     +     J>A     cosm^J    t+      5-B      sin  m    CO  t  (1) 

o  m  o  ^-    m  o 


where  f(t)  represents  the  time  series  being  analyzed,  m  is  the  integer 

specifying  the  particular  harmonic  being  calculated  and     LkJ     is  the 

o 

fundamental  angular  frequency,  which  is  related  to  the  time  duration  of 
the  series ,  T,  by 

The  power  present  in  the  mth  component  is  then 

P      =  1/2  (A    2  +  B    2)  (3) 

m         .       m  m 

and  the  magnitude  of  current  speed  in  that  component  is 

/      2  2,1/2 

C      =    A        +  B  (4) 

m  m  m 

Equivalently ,  Equation  (1)  may  be  expressed  as 

f(t)  =  C    cos(m   U3    t  -     Ot    )  (5) 

m  o  m 

where 

A     =  0  and      OC      =  sin"   Bm  (6) 

m 

Here      OC     is  the  phase  angle,  the  number  of  degrees  after  the  start  of 
m 

each  series  at  which  component  maximum  occurs. 

The  mean  velocity  and  the  phase  differences  between  the  diurnal 
component  and  the  transit  of  the  sun,  and  the  semi-diurnal  component 
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and  the  transit  of  the  moon  are  contained  in  Table  II,  along  with  the 
results  of  the  fourier  analysis.    The  total  power  present  in  all  the  cal- 
culable frequencies  is  also  tabulated  in  this  table  to  facilitate  a  com- 
parison of  this  value  with  the  amount  of  power  present  in  the  diurnal 
and  semi-diurnal  components.     The  conclusions  to  be  drawn  from  this 
table  are  the  following: 

(i)         The  amplitudes  of  the  diurnal  and  semi-diurnal  com- 
ponents are  comparable  with  the  mean  velocity, 
(ii)         The  power  in  the  diurnal  plus  semi-diurnal  components 

is  about  three-fourths  of  the  total  variable  power. 
(iii)         Diurnal  components  are  not  consistently  more  important 

than  the  semi-diurnal  and  vice-versa  . 
(iv)         The  variable  effects  are  generally  larger  in  the  E-W 

components  than  in  the  N-S  components,  while  for  the 
mean  velocities  the  reverse  is  true. 
(v)         The  diurnal  and  semi-diurnal  terms  in  both  the  N-S 

and  E-W  portions  of  the  current  have  a  fairly  consistent 
phase  relationship  with  the  meridional  transit  of  the 
sun  and  moon  respectively. 
The  phase  stability  of  the  fourier  components  may  be  seen  by 
comparing  the  phase  angles  relative  to  lunar  and  solar  transit  with  the 
respective  means. 

The  variations  from  the  mean  of  the  semi-diurnal  does  not  exceed 
+  2.5  hours  in  either  the  E-W  or  the  N-S  component  except  for 

67 


TA 

BLE  II. 

FOURIER  ANALYSIS  RESULTS 

se 

Mean  Vel. 
(cm/sec) 

Period 
(Hours) 

A 

B 

1  s 

O    CM 

Plj        — ' 
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P        ~ 

m 

J2  £ 
O   o 

P 

T 

q     Magnitude 
3        (cm/sec) 

q        Phase  Angle 
(degrees) 

Phase  Diff. 
(Hours) 

Crui 

Re  Sun 
Transit 

Re  Moon 
Transit 

24.8 

14 

-14 
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20.5 

-43 

-3.3 

I 

1.7 
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12.4 

-5 

4 

23 
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36 

+2.6 

24.8 
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1 

22 
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5.6 

51 

-2.3 

III 
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Cruise  III  in  the  E-W  and  Cruises  II  and  III  in  the  N-S.     Similarly,  the 
variation  from  the  mean  of  the  diurnal  phase  does  not  exceed  +  1.9  hours 
except  on  Cruises  V  and  VI  in  the  E-W  and  III  and  VI  in  the  N-S.     Only 
four  out  of  the  24  variations  from  the  mean  exceed  one-third  of  the  asso- 
ciated period.     Since  we  are  dealing  with  a  time  series  which  is  both 
non-stationary  and  short,  such  agreement  may  be  regarded  as  good. 
However,  the  fact  that  there  is  not  a  constant  relation  between  the 
amplitudes  of  these  two  components,  and  the  fact  that  the  phase  angles 
occasionally  vary  widely  from  the  mean  leads  to  either  of  two  inferences: 
(a)  the  currents  are  driven  indirectly  by  the  tides  and  are  considerably 
modified  by  some  unknown  oceanic  processes,  or  (b)  inertial  rotations 
are  important.     This  makes  one  unwilling  to  say  that  the  coupling  of  the 
tide  and  current  is  direct,  but  does  indicate  that  a  relationship  exists. 
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VII.     CONCLUSIONS 

Currents  in  Monterey  Bay  are  quite  variable  in  both  velocity  and 
direction.     Therefore,  one  cannot  describe  the  mean  flow  by  means  of  a 
short-term  measurement. 

All  six  cruises  yielded  mean  southerly  currents  .    The  average  of 
these  means  is  a  current  directed  163°T  at  a  speed  of  12.1  cm/sec 
(5.9  nm/day)  .     This  result  supports  Lammers'  analysis  that  the  California 
Current  sweeps  close  to  the  mouth  of  Monterey  Bay  (and  the  data  point) 
during  the  upwelling  period . 

It  has  been  established  that  there  are  substantial  diurnal  and 
semi-diurnal  components  which,  with  a  few  notable  exceptions,  have 
a  consistent  phase  relation  with  the  sun  and  moon.    There  obviously 
must  be  at  least  an  indirect  coupling  with  the  semi-diurnal  tide.     The 
diurnal  component  may  be  tidal  or  inertia  1  or  merely  related  to  the  passage 
of  the  sun . 
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APPENDIX  A 
COMPUTER  PROGRAMS 

1 .  Reduction  of  GEK  Data  to  Current  Vectors 

2.  Graphing  of  Ocean  Current  Data  versus  Winds  and  Tides 
3  .         Progressive  Vector  Diagram  of  Current  Vectors 

4.         Fourier  Analysis  of  Ocean  Current  Components 
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APPENDIX  B 
SAMPLE  DATA  FROM  CRUISE  I 
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CRUISE  I 
10  -  11  MAY  1972 


WIND  DATA 

TIME 

WIND  SPEED 

WIND  DIRECTION 

10  May  1972    1230 

8 

NW 

1415 

11 

WNW 

1600 

12 

WNW 

1715 

16 

NW 

1800 

20 

WNW 

1955 

14 

NW 

2100 

10 

NW 

2300 

8 

NW 

11  May  1972    2400 

8 

NW 

0200 

2 

VARIABLE 

0400 

4 

VARIABLE 

0800 

4 

VARIABLE 

0905 

6 

WNW 

1000 

7 

WNW 

1100 

8 

NW 

1200 

8 

WNW 

1300 

9 

WNW 
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CRUISE  I 
10  -  11  MAY  1972 


TIDE  DATA 

TIME 

TIDE  HEIGHT 

TIME 

TIDE  HEIGHT 

10  May  1972    1200 

-0.90                 11 

May 

1972 

2400 

-0.85 

1300 

-1.30 

0100 

-2.00 

1400 

-1.30 

0200 

-2.85 

1500 

-0.70 

0300 

-3.05 

1600 

+0.30 

0400 

-2.65 

1700 

+  1.35 

0500 

-1.80 

1800 

+2.40 

0600 

-0.70 

1900 

+3.15 

0700 

+0.40 

2000 

+3.40 

0800 

+  1.10 

2100 

+3.00 

0900 

+  1.45 

2200 

+2.00 

1000 

+  1.25 

2300 

+0.60 

1100 
1200 
1300 

+0.80 
+0.10 
-0.50 
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APPENDIX  C 
UNAVERAGED  CURRENT  DATA  PLOTTED  VERSUS  WINDS  AND  TIDES 
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